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The regions of hole concentrations < a; < 0.3 and temperatures 0.005|t| < T < 0.02|t| are studied 
in the t-J model of Cu-0 planes of perovskite high- Tc superconductors. For this purpose self-energy 
equations for hole and spin Green's functions are derived using Mori's projection operator technique 
and these equations are self-consistently solved. The calculated hole band transforms radically at 
X « 0.08. A narrow low-concentration band with minima near (±-|,±^) is converted to a band 
resembling the case of weak electron correlations, with the minimum at or (0,0). The hole 

Fermi surface is respectively changed from small ellipses at (±-|,±-|) to a large rhombus centered 
at (tt,?!) or (0,0). The decrease of the magnetic susceptibility at the antiferromagnetic wave vector 
and spin correlations with doping is determined by the growth of the frequency of spin excitations 
at this momentum. The shape of the frequency dependence of the susceptibility depends heavily on 
the hole damping and varies from a broad feature similar to that observed in La2-iSra;Cu04 to a 
pronounced maximum which resembles the normal-state resonance peak in YBa2Cu307-j;. 

PACS numbers: 71.10.Fd, 74.25. Ha, 74.25. Jb 



I. INTRODUCTION 

The two-dimensional t-J model was proposed by 
Anderson^ for the description of strong electron corre- 
lations in Cu-0 planes of perovskite high-T^ supercon- 
ductors. In Ref. |3 the similarity of the low-energy part 
of its spectrum with the spectrum of the realistic three- 
band Hubbard model was demonstrated. Nowadays the 
t-J model is one of the most frequently used models 
for the interpretation of experimental results in cuprates 
(for a review, see Ref. i. Different numerical and an- 
alytical methods were used for the investigation of the 
model. Among these methods are the exact diagonal- 
ization of small clustersf^ Monte Carlo simulations)^^ 
density matrix renormalization group calculations,^ spin- 
wave*'^ and mean-field slave-boson approximations.^*^ In 
spite of the considerable progress made towards the un- 
derstanding of the properties of the model, the basic is- 
sues of its behavior in going from light to heavy hole 
doping have not yet been completely resolved. In partic- 
ular, there is still no clear knowledge of how the narrow 
spin-polaron band inherent in light doping is converted to 
a wide band observed in photoemission of optimally and 
overdoped crystals^ The investigation of the variation 
of the magnetic susceptibility with the hole doping, tem- 
perature and the damping of excitations is also of great 
importance in view of a great body of data obtained in 
inelastic neutron scatteringpi^ 

Aiming at a description for the wide range of hole con- 
centrations in this paper we use the methoc^ based on 
Mori's projection operator techniqueii^ The method al- 
lows one to derive self-energy equations for Green's func- 
tions constructed from Hubbard operators without re- 
course to the intricate diagrammatic technique. These 
equations retain the rotation symmetry of spin compo- 
nents inherent in the Hamiltonian and do not imp ly any 
preset magnetic ordering. In contrast to Ref. [Ij, where 
this approach was implemented for the spin subsystem 



only, in this work we use it also for the hole subsystem. 



The obtained equations were self-consistently solved in 
a 20x20 lattice for the ranges of hole concentrations < 
X < 0.3 and temperatures 0.005|i| <T< 0.02|i| where 
t is the hopping constant of the t-J model. In cuprates 
this concentration range covers the regions from light to 
heavy doping. The ratios of the exchange J and hopping 
parameters of the model J/\t\ = 0.4 and J/\t\ = 0.2 were 
used. Such ratios correspond to hole-doped cupratesii^ 
For \t\ « 0.5 eV the boundaries of the above temperature 
range correspond approximately to 30 K and 120 K. 



Obtained results indicate that the hole band trans- 
forms radically at x ~ 0.08. A narrow low-concentration 
band with minima near (±^, ±-|) is converted to a band 
resembling in its shape the case of weak electron correla- 
tions, with the minimum at (7r,7r) or at (0,0) in depen- 
dence on the sign of t. The hole Fermi surface is respec- 
tively changed from small elliptical pockets at (±|, ±|) 
to a large rhombus centered at (tt, tt) or (0,0). With in- 
creasing X the maximum in the imaginary part of the 
magnetic susceptibility x" ^-t the antiferromagnetic wave 
vector Q = (tt, tt) loses its intensity and shifts to higher 
frequencies. This behavior is mainly connected with the 
growth of the frequency of spin excitations at Q and is 
consistent with experimental data. The shape of the fre- 
quency dependence of x"(Q) depends heavily on the hole 
damping and varies from a broad feature similar to that 
observed''^- in La2-a;Sr2;Cu04 to a pronounced maximum 
which is analogous to the normal-state resonance peakii 
in YBa2Cu307_y. 
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II. MAIN FORMULAS 

The Hamiltonian of the two-dimensional t-J model 
reads 



H 



nmcT 



t 



2 ^ 

nm 



V''n''m 



(1) 

where a^a — Incr) (nO| is the hole annihilation operator, n 
and m label sites of the square lattice, ct = ±1 is the spin 
projection, Inu) and |nO) are site states corresponding to 
the absence and presence of a hole on the site. These 
states may be considered as linear combinations of the 
products of the ?>dx^~y^ copper and 2p^ oxygen orbitals 
of the extended Hubbard modelii^ We take into account 
nearest neighbor interactions only, = t <5n,m+a 
and — J X^a i^n.m+a where the four vectors a con- 
nect nearest neighbor sites. The spin-i operators can be 
written as = ^ tT|ncr)(n(T| and — \r\a){n,—a\. 

Properties of the model are determined from the hole 
and spin retarded Green's functions 



G(kt) = -,:0(t)({ak.(<),aL}), 



(2) 



D(ki) = -*0(i)([4W,5-k]>, 



where Oko- and are the Fourier transforms of the re- 
spective site operators, operator time dependencies and 
averaging are defined with the Hamiltonian H = H — 



with the chemical potential ii. As men- 



tioned, to obtain self-energy equations for these func- 
tions we used Mori's projection operator technique^^ii^ 
In this approach the Fourier transform of Green's func- 
tion {{Ao\AI)) is represented by the continued fraction 



Poi4)) = 



1^0-41 



U} - En - 



Vn 



(3) 



The elements of the fraction Ei and Vi are determined 
from the recursive procedure 



[An,H] — EnAn + ^n+l + T4_iA„_i, 

E„^\[A„,H]-Al\\A„-Ai\-\ 

Vn-i = \ An ■ Ai\ |A„_i • aI_-^\-\ 
V-i=0, n = 0,1,2,... 



(4) 



The operators Ai constructed in the course of this proce- 
dure form an orthogonal set, \Ai • oc Sij. In Eqs. Q 

and Q the definition of the inner product \Ai ■ Aj\ 
depends on the type of the considered Green's func- 
tion. For example, for functions Q these are {{Ai,Aj}) 
and ([Ai,^]]), respectively. The method described by 
Eqs. (O and can be straightforwardly generalized 



to the case of many-component operators which is nec- 
essary, for example, to consider Green's functions for 
Nambu spinors in the superconducting state>i^ 

The residual term of fraction ||2J) is the Fourier trans- 
form of the quantity 



T=|A„,.At||A„_i.4_J-i 



(5) 



where the time evolution of the operator An is deter- 
mined by the equation 



dt 



Ant = Y[il - Pk)[Ant,nl An,t=a^An (6) 



fe=0 



with the projection operators P„ defined as PnQ = 
\Q ■ Al\\An ■ Al\-^ An. The residual term T IS a many- 
particle Green's function which can be estimated by the 
decoupling. Following Ref. the decoupling procedure 
may be somewhat improved by introducing the vertex 
correction a which is determined from the constraint of 
zero site magnetization 

{s'i.) = \{l-^)-{s;.'st')=Q. (7) 

The spin Green's function obtained in this way readsi^ 

4(7k-l)(JCi+tFi) 



D(ku;) 
Imn(kw) 



w2 - tjn(kcj) - u^l ' 

<dTTej^{l-x) 
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where 7k = ^[cos(fcj;) -t- cos(fcj^)], Ci = (sit^^'^n+a) ^^d 
^1 — (anOn+a) are the spin and hole correlations on 
neighboring sites which can be derived from the respec- 
tive Green's functions, 
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(4'^m> = T^E^^''^"'™^ / dco[l + nB{u;)]B{ku;), 

(9) 

daj7T,i?(aj)A(kti;), 



k 



= 16j2a|Ci|(l-7k)(A + l + 7k), 



(10) 



is the square of the frequency of spin excitations with the 
parameter A describing the spin gap near Q, A'^ is the 
number of sites, A(ka;) — — 7r~^ImG(k(jj) and B{kLij) — 
— 7r~^ImD(kti;) are the hole and spin spectral functions, 
npii^) = [exp(a;/T) + l]-i, nsiy^) = [exp(c<j/T) - 1]"! 
with the temperature T . The real part of the self-energy 
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n(kci;) can be calculated from its imaginary part and the 
Kramers-Kronig relation. However, in the calculations 
Re n was neglected, since the gap parameter A and along 
with it the magnon frequencies near Q were determined 
from the constraint ((TJ. This momentum region plays the 
main role in the considered spectral and magnetic prop- 
erties. As follows from earlier calculationsji^ the vertex 
parameter a depends mainly on x and can be approxi- 
mated as a = 1.802 - 0.802 tanh(10a;). 

Two steps of the procedure (|3J| , Q can be also carried 
out for the hole Green's function. We find 

({aka,aL}) = + ^) 

Eo = {4t(f> + QtCi(t>-^ - 3 JFi?!)-i)7k 



(11) 



^x{l - x) 



4Ci 



+ 24<2Fi7k - 4t(3Ci0-i + 2cj))iEo + fihi. 
El w 

where x = N-^ J2u IZo dujnF{Lo)A{kuj) and C^^ ^ 2N-^ 

Sk' 7k-k;'('5k"*-k')- Continued fraction ||2Jl is ter- 

minated at this stage, the hole spectrum is approximated 
by two poles which form two bands, 



G(kcj) = 



- ^) / 1 

1 



uj - £k,2 + l-J. + iri 



(12) 



where 



Eo = {4t4> + mCi4)-^ - 3JFi0-i)7k, 



(13) 



Ek, 



Eo/2±{El/^ + Vofl\ 



and 77 is a damping of the hole states. In the framework 
of the t-J model this damping is connected with the hole- 
magnon scattering. However, in this work the damping 
is considered as a free parameter to take into account 
other possible damping processes. Such approach is mo- 
tivated by the fact that the damping of spin excitations 
(|H1) depends heavily on the value of 77 and to investigate 
different shapes of the frequency dependence of the sus- 
ceptibility 77 is considered as a variable parameter. 

Thus, in the present work the calculation of the ele- 
ments of the continued fraction Ei and Vi for the spin 
Green's function was carried out up to the second order: 
Eq , Vo and Ei were calculated. Then the residual term of 
the continued fraction (the part of the fraction containing 
the third and higher order terms) which is represented by 
the many-particle Green's function (O is approximated 



by the decoupling. The same elements of the continued 
fraction were calculated for the hole Green's function [see 
Eq. l(TT|l ]. In this case the residual part of the continued 
fraction was omitted and the artificial damping 77 was 
added. The main argument for the above truncations of 
the continued fractions is based on the results obtained 
previously. In Ref. the analogous truncation in the 
equations of motion with the decoupling improved by a 
vertex correction was successfully used for the descrip- 
tion of spin excitations in the Heisenberg model. For 
the lightly doped t-J model this approximation was used 
in Ref. 113 and the results obtained were in good agree- 
ment with the exact diagonalization of small clusters. 
The two-pole approximation for the hole Green's function 
was shown to give a good description of the spin-polaron 
band for moderate \t\/ J in the spin-wave approach.* In 
a more general perspective the used method is a version 
of the Lanczos algorithm. It is known from numerical 
methods that at a successful choice of the starting vector 
this algorithm gives good approximations for outermost 
eigenvalues and eigenvectors of a matrix already with few 
steps. 

Equations (f7jl- (fT^ form a closed set which can be 
solved by iteration for given values of the chemical po- 
tential, temperature, hole damping and the sign of t. 



III. EVOLUTION OF THE HOLE SPECTRUM 

Results of our calculations for the lowest of the hole 
bands which crosses the Fermi level are shown in Figs. 
and El for different values of the hole concentrations and 
the sign of t. Notice that for low hole concentrations the 
shape of the band is close to the shape of the spin-polaron 
band obtained in the spin- wave approximation3*2, and in 
the approach of Ref. 1131 However, for the parameters 
chosen the bandwidth is half as much again the value 
found in the two mentioned approaches. This difference 
is a consequence of our two-pole approximation H12|l for 
the hole Green's function. For small x the shape of the 
band depends only weakly on the sign of t. With increas- 
ing X the band changes rapidly and already for moderate 
concentrations its shapes for opposite signs of t differ es- 
sentially and resemble shapes of weakly correlated bands. 
As this takes place, the bandwidth grows and for x w 0.3 
it is three times as much as the bandwidth for light dop- 
ing. 

The decrease of the temperature to T = 0.005|i|, of 
the hole damping to 77 = 0.015|t| or the increase of the 
exchange constant to J = 0.4|t| do not change qualita- 
tively the evolution of the hole band with doping. For 
all these parameters the transformation of the hole band 
from the shape typical for strong correlations to that re- 
sembling the weakly correlated case occurs at cc ~ 0.08. 
As this takes place, the respective Fermi surface, shown 
in the base planes in Figs.nandlJl changes its shape from 
small ellipses at (±-|, ±^) for small a; to a large rhombus 
centered at (tt, tt) or (0, 0) for large x. The van Hove sin- 
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77 = O.OSi and x = 0.033, 0.076, 0.14, and 0.18 (from top to 

bottom). Horizontal planes arc the Fermi levels. The respec- FIG. 2: The hole dispersion for t < 0, J/\t\ = 0.2, T = 0.02|tl, 
tive Fermi surfaces are shown in the base planes. n = Q.05\t\ and x = 0.033, 0.076, 0.13, and 0.2 (from top to 

bottom). 
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gularity of the hole band located at (0, tt) remains near 
the Fermi level in the wide range of hole concentrations 
X < 0.2. 

The reason for the above transformation of the band 
shape is in the following. For small hole concentrations 
the hole dispersion is determined by the magnetic order- 
ing which is close to the long-range antiferromagnetic or- 
der. In this case the elementary cell in the direct space is 
doubled due to the opposite alignment of spins on neigh- 
boring sites. Hence the Brillouin zone becomes half the 
value of the paramagnetic zone and the points (0, 0) and 
(tt, tt) become equivalent. In particular, this means that 
the hole dispersion is approximately invariant with re- 
spect to the substitution k — > k -f Q, as seen in the upper 
parts of Fig.nand|21 Mathematically this corresponds to 
the case when |£'o | <C Vb in £k j- in Eq. ifT^ . With increas- 
ing X the spin correlation length ^ is decreased. Whenii^ 
f = iflA"^/^ becomes comparable to a few lattice spac- 
ings a the above arguments cease to work. In this case 
the hole dispersion is determined by Eq which becomes 
larger than Vb and the band acquires the shape which is 
similar to the case of weak correlations. However, even 
for X = 0.25 the bandwidth is much smaller than 8|i|, 
the bandwidth in the latter case (see the lower parts of 
Fig. n and 121). Thus, even in the overdoped case the cor- 
rections due to electron correlations are essential. The 
above transformation is not connected with some sharp 
transition in the magnetic subsystem. The decrease of 
the spin correlation length with x is continuous. 

The analogous transformation of the hole band is ex- 
pected also in the more complicated and flexible t-t'-t"-J 
model which is frequently used for the interpretation of 
experimental results in cuprates^*ii The respective value 
of the hole concentration is supposed to be close to 0.08. 
This assumption is based on the fact that the transforma- 
tion occurs when the correlation length approaches a few 
lattice spacings. The correlation length derived from neu- 
tron scattering experimentslS in cuprates acquires this 
value near the mentioned concentration. 

In the considered two-pole approximation only one hole 
band crosses the Fermi level. In contrast to this the 
calculations based on the spin-wave approximation^ and 
on the approach of Ref. [ij, which are supposed to be 
more accurate for low hole concentrations, show that 
for X > 0.05 there are two bands crossing the Fermi 
level. One of these bands is the spin-polaron band which 
corresponds to a pronounced peak in the hole spectral 
function and has the dispersion similar to that shown 
in upper parts of Figs. and |21 The second band has 
a much weaker and broader peak in A{\iuj) and there- 
fore its contribution to the lower band of the two-pole 
approximation can be neglected. Thus, the lower hole 
band which cross the Fermi level corresponds to the fea- 
ture with the highest peak intensity in the hole spectrum 
- the spin-polaron band. The upper band corresponds to 
"everything else" in the hole spectrum. From previous 
considerations (see Refs. I^Sllj and references therein) 
it is known that in A(kw) for low x several weaker and 



broader maxima are located above the spin-polaron peak. 
These maxima merge into one broad maximum at larger 
concentrations. The upper band of the two-pole approxi- 
mation gives a rough description for these peculiarities of 
the hole spectrum. Since states near the Fermi level are 
of main interest and this approximation gives a satisfac- 
tory description for these states, it is used in the present 
work. 

It is worth noting that the used two-pole approxima- 
tion for the hole Green's function gives satisfactory re- 
sults not only in the limit of small x but also in the limit 
X 1. Indeed, in this limit Vq ^ 0, Eq k 4t7k — /i 
and GCkuj) — [uj — Atji^ + fj,)~^. However, the two-pole 
approximation becomes inapplicable in the intermediate 
region for a; > 0.3 when the chemical potential falls into 
the gap between two bands. 



IV. MAGNETIC SUSCEPTIBILITY 

The imaginary part of the magnetic susceptibility is 
given by the equation 

X"(kw) = 47r^|B(kw) 

^ 16^^1(1 - 7k)(JCi +^j^iMniH 

{lu^ -LuRell-ujl)^ + {ujlran)^' ^ ' 

where fiB is the Bohr magneton. For the considered hole 
concentrations the hole contribution to the susceptibility 
is much smaller than the spin contribution.^ Up to x w 
0.13 the momentum dependence of x" is strongly peaked 
at Q (the use of a comparatively small 20 x 20 lattice does 
not allow us to describe the incommensurability of the 
magnetic response - the low-density inverse space and the 
frequency independent damping 77 produce a minimum in 
the spin-excitation damping ImH(k) at k = Q which is 
too shallow to give rise to incommensurabilit}ii2iS&) . The 
calculated frequency dependencies of the susceptibility at 
this momentum are shown in Fig.|2| In Lai.8gSro.i4Cu04 
the susceptibility is peaked at incommensurate momenta 
k = (tt ± 27r^, tt), (7r,7r ± 27r^). The experimental sus- 
ceptibility measuredi^ at one of these momenta is also 
shown in Fig. To demonstrate the similarity of the 
experimental and calculated dependencies the former was 
increased by approximately 2.5 times (the difference in 
the peak amplitudes is apparently connected with the 
splitting of the commensurate peak into 4 incommensu- 
rate maxima). As mentioned above, for the comparison 
with experiment we set t = 0.5 eV, thus the tempera- 
tures T — 0.005t and T — 0.02t in the figure correspond 
to 29 K and 116 K, respectively. 

As seen from Fig|2K, the calculated frequency depen- 
dencies are close to those observed in Lai.86Sro.i4Cu04 
for both temperatures. For T = O.OOSt the susceptibil- 
ity has a broad maximum at a; « 7 meV. This maxi- 
mum is not connected with the resonance denominator in 
Eq. H14|) . because for the parameters chosen spin excita- 
tions near Q are overdamped - their dampings are larger 
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FIG. 3: The frequency dependence of the imaginary part of 
the susceptibility, (a) The sohd hne corresponds to J = 0.2t, 
t > 0, T = O.OOSt, T] = 0.05t and x ^ 0.1. The dashed hne 
are for T = Q.02t and the same other parameters. Symbols 
are experimental results— in Lai.86Sro.i4Cu04 at T = 35 K 
(filled circles) and T = 80 K (open circles). Vertical bars show 
experimental errors. The calculated susceptibility is given for 
k = Q, the experimental data are for the wave vector of the 
incommensurate peak, (b) The solid line is the same as in 
part (a). The dashed line is calculated for 77 = O.OlSf, other 
parameters are the same as for the solid line. 



than their frequencies. A completely different situation 
occurs with the decreased hole damping. The dampings 
of spin excitations decrease also, the excitations cease to 
be overdamped and the shape of the susceptibility is de- 
termined by the resonance denominator in Eq. 114|l . In 
this case x" has a pronounced maximum at w w Wk, as 
seen in Fig. ^p. The shape of the dashed curve in this 
figure is similar to the resonance peak observed in the 
normal state of the underdoped YBa2Cu307_yi^^ Thus, 
we suppose that the dissimilarity of the frequency de- 
pendencies of the susceptibility in this crystal and in 
La2_a;Sra;Cu04 is connected with different values of the 
damping of spin excitations. One of the possible reasons 
for this difference is the diverse values of the hole damp- 
ing. 

The evolution of x"{Q^) with the increase of the hole 
concentration is shown in Fig. 01 For the damping chosen 
all 3 curves have pronounced maxima, however also for a 
larger damping, when spin excitations are overdamped. 
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FIG. 4: The susceptibility at the antiferromagnetic wave vec- 
tor for the 3 hole concentrations indicated near the curves. 
Other parameters are J = 0.2t, T = 0.02i, rj = 0.02t and 
t > 0. 



the behavior of the susceptibility is similar to that shown 
in the figure. With growing x the frequency of the maxi- 
mum increases and its amplitude decreases. This behav- 
ior points to the increase of the spin excitation frequency 
ujQ with the hole concentration. It can be showniS that 
the value of this frequency is directly connected with the 
magnetic correlation length ^. From the obtained re- 
sults it follows that for low hole concentrations and tem- 
peratures ^ ~ ax~^/^ where a is the intersite distance. 
An analogous relation has been derived from experimen- 
tal data in La2_a:Sr2,Cu04iSi The peak amplitudes of 
the susceptibility decreases rapidly with increasing x for 
X ^ 0.11 and then flattens out. The rapid decrease can be 
related to the transformation of the hole band discussed 
in the previous section. The qualitative behavior of the 
peak amplitude is not changed with the increase of the 
exchange constant to J = OAt, with the decrease of the 
temperature to T = O.OOSt or with the change of the sign 
of t. The qualitatively similar behavior of the suscepti- 
bility is observed'^^ in the normal-state YBa2Cu307_y. 
However, in experiment the decrease of the peak inten- 
sity is somewhat slower for small x and the maximum 
of the susceptibility is smeared out for x > 0.16. These 
features depend strongly on the damping of spin excita- 
tions which is described only crudely in the used two-pole 
approximation. 

The decrease of the susceptibility at the antiferromag- 
netic wave vector with doping is reflected in the weak- 
ening of spin correlations Ci — {s^^Sq^). These correla- 
tions are shown in Fig. [S] where the components of the 
vector 1 are designated as n and m. In this figure the cor- 
relations obtained^ by exact diagonalization in a 4x4 lat- 
tice for J = 0.4i and T = are also shown. The difference 
between the two sets of data is connected with the influ- 
ence of finite-size effects in the latter data and with the 
known issue of the method used in our calculations which 
somewhat underestimates the correlations.^- '^^ The cor- 
relations change only slightly with the change of the sign 
of t and with the variation of J and T in the considered 
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0.00 0.05 0.10 0.15 0.20 0.25 
X 

FIG. 5: The spin correlations Ci = {s^^s'^^), 1 = {n,m) 
obtained in our calculations for J = OAt, T = 0.005f, 
77 — 0.05t and t > (solid lines) and in exact-diagonalization 
calculations'' in a 4x4 lattice for J — OAt and T — (sym- 
bols, dashed lines are drawn as a guide for the eye). 

ranges. In spite of the large sensitivity of the spin suscep- 
tibility on r], the spin correlations depend only weakly on 
the hole damping due to the integration over frequencies 
in the formula for the spin correlations ©. The change 
of the hole damping from 77 = 0.05t to 0.015t leads to the 
variation in Ci by several percent. 

V. SUMMARY 

In this paper the evolution of the hole and spin- 
excitation spectra of the two-dimensional t-J model was 
studied in the range of hole concentrations < a; ^ 0.3 



which spans the regions from light to heavy doping. The 
variation of the spectra with temperature, with the sign 
of the hopping parameter and with the excitation damp- 
ing was also considered. For this purpose the self-energy 
equations were derived employing the projection opera- 
tor technique and these equations were self-consistently 
solved. The hole band was found to transform radically 
at a; « 0.08. A narrow low-concentration band with 
minima near (±^, ±-|) is converted to a wider band re- 
sembling in shape the case of weak electron correlations, 
with the minimum at (tt, tt) or at (0, 0) in dependence 
on the sign of t. The hole Fermi surface is respectively 
changed from small elliptical pockets at (±|-,±-|) to a 
large rhombus centered at (7r,7r) or (0,0). The frequency 
dependence of the imaginary part of the magnetic sus- 
ceptibility x" ^-t the antiferromagnetic wave vector de- 
pends heavily on the damping of spin excitations and 
varies from a broad feature similar to that observedji^ in 
La2_2;Sr2;Cu04 to a pronounced maximum which resem- 
bles the normal-state resonance peak in YBa2Cu307_y 
One of the possible reasons for the variation of the spin- 
excitation damping is the change of the hole damping. 
With increasing doping the maximum in the suscepti- 
bility loses its intensity and shifts to higher frequencies. 
The similar behavior is observed in cuprates and is con- 
nected with the growth of the spin-excitation frequency 
at the antiferromagnetic wave vector, which reflects the 
decrease of the magnetic correlation length with doping. 
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